Topological Frustration in Graphene Nanoflakes: 
Magnetic Order and Spin Logic Devices 
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Magnetic order in graphene-related structures can arise from size effects or from topological frus- 
tration. We introduce a rigorous classification scheme for the types of finite graphene structures 
(nano-flakes) which lead to large net spin or to antiferromagnetic coupling between groups of elec- 
tron spins. Based on this scheme, we propose specific examples of structures that can serve as 
the fundamental (NOR and NAND) logic gates for the design of high-density ultra-fast spintronic 
devices. We demonstrate, using ab initio electronic structure calculations, that these gates can in 
principle operate at room temperature with very low and correctable error rates. 



Magnetism in solids is typically associated with the 
presence of transition metal elements with d electrons. 
Experimental evidence has recently shown that mag- 
netism can also arise in solids composed of elements with 
only sp electrons, as, for example, in bulk proton irra- 
diated graphite [H]. Theoretical studies addressing the 
issue of magnetism in carbon-based materials have fo- 
cused onpoint defects 0, Q and reduced dimension- 
ality [i, [1, 0, S, [i, These works have shown that 
magnetism can arise in various situations, for example, 
the antiferromagnetic order across the edges of zigzag- 
edged graphene nanoribbons [5] and the large net spin in 
zigzag-edged triangular graphene nanoflakes 0, B Q ■ 
the other hand, magnetism from sp electrons is not new 
in organic chemistry: high-spin states have been observed 
in synthesized organic molecules, for example, conjugated 
polyradicals [lHand triangulene derivatives where 
all conjugated tt electrons can not be paired simultane- 
ously. At present, it is not yet clear if the magnetic prop- 
erties of irradiated bulk graphite and those of high-spin 
organic molecules have the same origin. Understand- 
ing of the magnetic order and its origin in various fi- 
nite graphene structures is of fundamental importance 
as well as of practical interest for application in spin- 
tronics where carbon based materials have recently 
demonstrated their potential 14 1. 

In this letter, we show that a distinct origin of mag- 
netism in finite graphene structures is topological frus- 
tration of TT bonds, which is a generalization of the sim- 
ple counting rule that governs magnetic order in organic 
molecules. We use the notion of this topological frustra- 
tion to derive a rigorous classification scheme for arbitrar- 
ily shaped graphene nanoflakes, depending on whether 
only one or both sublattices of the graphitic structure 
are frustrated. From this classification scheme we iden- 
tify which nanoscale structures can give rise to a strong 
antifcrromangetic (AF) coupling and propose a specific 
example of a structure which can serve as the fundamen- 
tal (NOR and NAND) spin logic gate. Finally, we employ 



first-principles electronic structure calculations to show 
that this type of gate can operate at room temperature, 
an important prerequisite for the design of realistic, high- 
density ultra-fast spintronic devices based on graphene. 

A graphene nanoflake (GNF) is an arbitrarily shaped 
finite graphene fragment consisting of hexagonal rings 
and bounded by a single (non-self-intersecting) topolog- 
ical circuit, where all in-plane dangling cr-bonds at the 
edge are assumed passivated. We start the general clas- 
sification of GNFs with the widely used Pz band Hub- 
bard model, where magnetic correlations are described 
through the hamiltonian 
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In the first term, the tight-binding part, the operators 
Cia (cjg.) annihilate (create) an electron at site i with 
spin a =t,J, and t is the hopping integral between the 
nearest neighbor sites i and j. The honeycomb lattice of 
graphene is bipartite, that is, any pair of bonded nearest 
neighbor carbon atoms consists of one atom from each 
of the two interpenetrating sublattices, commonly de- 
noted as A and B (Fig. 1). The second term of the 
hamiltonian describes electron-electron interactions via 
the on-site Coulomb repulsion U with = c\^Cia the 
number operator. This interaction may trigger an in- 
stability in the low-energy electronic states and produce 
spin polarized states to minimize the total energy. Zero- 
energy (non-bonding) states in a half-filled tt sub-band 
are particularly prone to become polarized for U > Q. 

The occurrence of zero-energy eigenstates in the tight- 
binding hamiltonian for a GNF can be accounted for by 
a theorem on hexagonal graphs [l5| . The number of such 
states, called "nullity" (77), is determined by the topology 
of the GNF according to the equation rj ~ a — (3, where 
a and (3 are the maximum numbers of non- adjacent ver- 
tices and edges respectively. The latter is also called the 
maximum matching of the GNF graph and satisfies the 
relations (i — 6 — v, where 9 and v are the numbers of 
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FIG. 1: (Color online) (a) and (b) Class I GNFs, where nul- 
lity is equal to sublattice imbalance, i] = |A'a — A'bI- (c) 
Class II bowtie-shaped GNF with zero sublattice imbalance 
but a nullity of two, specifically, A'a = A'b = 19, a = 20. 
The definition of a requires switching of sublattice across the 
dashed line. All colored sites correspond to a maximum set 
of non-adjacent sites. 

positive and negative eigenstates respectively. The sum 
a + /3 equals to N, the total number of carbon atoms 
in the GNF. When 77 = 0, /3 = N/2 and all carbon 
atoms can be connected by a set of non-adjacent pair- 
wise bonds, which is referred to as "perfect matching" 
indicating a perfect pairing of all orbitals. Otherwise, 
rj — a~f3^N — 2f3>0, which is the number of sites 
that are left out by the best possible matching. The in- 
ability to simultaneously pair all Pz orbitals is entirely 
attributed to the topology of the GNF, therefore can be 
called topological frustration. 

In general, GNFs can be classified into two classes ac- 
cording to whether one or both sublattices are topolog- 
ically frustrated. In class I, at most one of the sub- 
lattices is frustrated, which is characterized by /3 = 
minjiVA, Nb}-, that is, the maximum matching covers all 
sites of at least one sublattice. This class includes all 
highly symmetric forms of GNFs as we proved previously 
0. The nullity for this class is simply rj = \Na — Nb\, 
as illustrated in Fig. [ija) andfTJb). Balanced sublattices 
{Na = Nb) means zero nullity. In class II, both sublat- 
tices are frustrated, characterized by /3 < minjA^A, A'a} 
(note (3 < minjiVA, A^b}), which means r] > \Na — A^b|- 
As a result, the nullity can be finite even for GNFs with 
balanced sublattices. An example is shown in Fig. [Ijc), 
where the nullity is 77 = 2 even though Na = Nb- 

Although graph theory is all that is required to pre- 
dict the number of singly-occupied orbitals, it is not 
clear how the electron spins in these orbitals are aligned. 
Complementary information comes from the Lieb the- 
orem [3| which determines the total spin but not the 
number of singly-occupied orbitals. The Lieb theo- 
rem was proved for any even-numbered bipartite sys- 
tem, where the ground state has a total magnetic mo- 
ment S ^ \Na~ Nb\/2. Therefore, for class I GNFs, 
S = \Na ~ Nb\/2 = 77/2, that is, all spins in singly- 
occupied orbitals align parallel to each other, consistent 
with Hund's rule. This was confirmed experimentally [12] 
and by first-principles calculations 0, Q. For class II 



GNFs, 5" = |7Va - iVB|/2 < 77/2, indicating the existence 
of antiferromagnetic (AF) order, and Hund's rule breaks 
down. Specifically, for the bowtie shaped GNF shown 
in Fig. [TJc), the magnetic moments of the left and right 
triangle must be AF coupled to satisfy the requirement 
5 = 0. This is also consistent with the fact that magnetic 
moments are localized in the two sublattices of graphene 
favoring AF coupling [1, [l3| . 

Besides topological frustration, AF coupling can also 
be induced by the polarization of the low-energy states 
that approach the Fermi level as the system size in- 
creases. This is a distinctly different origin of magnetism, 
for two reasons: first, it can not give rise to net spin; sec- 
ond, the energy of the non-interacting eigenstate is not 
strictly at the Fermi level, except in infinite systems, and 
magnetic order appears only if the interaction energy U 
is above a positive threshold or, equivalently, if the sys- 
tem is above a critical size. Examples of this mechanism 
are graphene nanoribbons [l8l | and hexagonal graphene 
nanoflakes 

We next focus on the magnetic coupling induced by 
topological frustration in class II GNFs and use the 
bowtie structure of Fig. [T]^c) as the simplest represen- 
tatives. Such structures have a low-spin ground state, 
involving spins spatially segregated and AF coupled. 
This open-shell, low-spin feature is not only fundamen- 
tally interesting but also may enable practically acces- 
sible logic operations. For instance, the simple bowtie 
structure is a natural NOT gate because flipping the in- 
put spins on one side of the bowtie requires the output 
spins on the other side to flip as well since the spins 
on the two sides must point in opposite directions as 
long as the AF order is the ground state. Practically, 
various means may be used to flip the input spin, in- 
cluding polarized light, local magnetic fields or direct in- 
jection of polarized electrons through magnetic materi- 
als. Among those, spin injection [l9!| appears the most 
promising, especially when considering the natural inte- 
gration of GNFs with graphene nanoribbons which have 
been predicted to exhibit rich spintronic properties in- 
cluding half metallicity 0, Q. The barrier for flipping 
the output spin is expected to be extremely low (^ ^bT) 
due to the weak spin-orbit coupling in carbon materi- 
als , which is a prerequisite for operation with low en- 
ergy consumption. In order to achieve both ultra-fast 
switching and robust operation at reasonable tempera- 
tures, it is also important to have the magnetic cou- 
pling 2 J = EpM — Eaf, the energy difference between 
the ferromagnetic and antiferromagnetic conflgurations, 
be greater than 18 meV, the minimum energy dissipa- 
tion [231 A;rBln2 evaluated at 300 K. In principle, pi- 
cosecond flipping of electron spins can be achieved with 
an energy splitting J > 7r?i/10^^^ sec ~ 2 meV. 

In order to establish whether or not the above condi- 
tions can be met, we turn to flrst-principles electronic 
structure methods to investigate the energetics of mag- 
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FIG. 2: (Color online) (a) The spectrum of singly occupied 
states of a bowtie-shaped GNF with spin up (lu and 2u) 
and down (Id and 2d), (b) Isodensity surface of the total 
spin distribution showing opposite spins localized at opposite 
sides, (c) Wavefunctions of the four singly occupied states. 
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FIG. 3: (Color online) (a) Dependence of the spin-coupling 
2J on the GNF bowtie geometry for (b) symmetric and (c) 
asymmetric configurations of various triangular sizes n and 
junction widths m. The dashed horizontal line marks the 
minimum coupling required for room temperature operation. 



netic coupling in detail. Calculations were performed us- 
ing spin-polarized density functional theory (DFT) as im- 
plemented in the SIESTA code [2lj. The generalized gra- 
dient approximation (GGA) exchange-correlation func- 
tional ^2^] was employed together with a double-^ plus 
polarization basis set, norm-conserving pseudopotentials 
[23'| and a mesh cutoff of 200 Ry. The different spin con- 
figurations were obtained by means of providing appro- 
priate initial guess electron spin densities 2J, |25| . The 
electronic structure of a representative bowtie GNF is 
shown in Fig. [21 The nullity of the corresponding graph, 
?7 = 4, manifests itself as four singly-occupied orbitals. 
Within the spin-polarized formalism, when spin-spatial 
symmetry breaking is allowed these four states are split, 
as shown in Fig. [Ifa). The wave functions clearly illus- 
trate the open shell singlet nature of the system: elec- 
trons are AF coupled but not paired in the sense that 
they are spatially segregated, in agreement with previous 
calculations J2d\ on non-kekulean molecules. The spin 
coupling here is 2J — 45 meV, well beyond the above 
mentioned thermodynamic threshold. 

The coupling strength in a bowtie GNF can be fur- 
ther engineered by optimizing its geometry. To explore 
possibilities, we calculated bowtie nanoflakes of different 
triangle size n and junction width m (both measured in 
units of the graphene lattice constant, oq = 0.25 nm) 
for both symmetric and asymmetric geometries, shown 
in Fig. [21 In either case, there are n — to — 1 non-bonding 
states on each triangularly shaped side of the bowtie 
GNFs. The calculated coupling magnitude 2J converges 
quickly with increasing n after reaching a maximum value 
Fig. [3^a). The leveling off is attributed to the size effect 
which diminishes the minimum energy splitting. Inter- 
estingly, for a specific asymmetric to = 2 configuration, 
the interaction strengths are as large as 2J = 180 meV. 



In comparison, the AF cou plin g of quantum dots [27 1 
and transition metal atoms [28| suffer from weak maxi- 
mum coupling strength, about 1 meV and 6 meV respec- 
tively, limiting their operation to very low temperatures. 
With a coupling of 180 meV, a GNF-based spin gate 
may operate at room temperature with an error rate of 
p — e^^'^/'^^B — 0.001, which can be handled by error 
correction schemes. 

To further illustrate the concept of using Class II 
GNFs for spin logic processing, we explore a tri-bowtie 
GNF in which the central triangle is connected to three 
other surrounding triangles through its vertices, as shown 
in Fig. [HJa). From graph theory arguments, we ex- 
pect the number of unpaired spins in the central re- 
gion D to be nu — toad — wbd — "t-cd — 1 where mxY 
{X, Y = A,B, C, D) are the widths of the junctions; the 
peripheral triangles accommodate ua — toad^I unpaired 
spins. The unpaired spins in the central region now de- 
pend on the competition of the spins of all the three 
peripheral regions: they tend to be AF coupled to the 
majority spins of the peripheral regions so that the total 
energy is lowered. We can therefore assign the total spins 
of two of those peripheral regions as operands A and B 
and the third as a programming bit C; the spin in the 
central region D is the output. With the spin- up state 
representing 1 and the spin-down state 0, the above logic 
is written as 



D = {Ar\B)u{Br\C)u{CnA) 
= {AnB)u{{AuB)nC) 



If C is 1, the logic reduces to D 
if C is 0, 13 



AUB, a, NOR gate; 



An B, a, NAND gate; these are the two 
fundamental gates in Boolean logic, from either one of 
which all other gates can be constructed. Therefore, the 
above scheme in principle provides an efficient design, 
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design are not beyond the reach of modern nano-scale 
fabrication methods. 
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FIG. 4: (Color online) (a) Reconfigurable spin logic NOR 
and NAND gate based on of a tri-bowtie GNF structure with 
UA = na = nc = 4, WD = 6, m = 1 [A, B and D are two in- 
puts and one output, respectively, and C is the programming 
bit), (b) A scheme of the localized spins and the couplings 
(2JxY = 34 meV). (c) Two distinct spin configurations corre- 
sponding to 1110 and 0110 for the ABCD spins, respectively, 
(d) The truth table of the programmable logic gate and the to- 
tal energy Etot of the operation configuration; D' and E[^^ are 
the error output and the corresponding energy {E[at > -Btot). 



with the spin degree of freedom employed to satisfy in 
the classic regime all digital logic operations. 

Undoubtedly, various engineering issues will have to be 
addressed before the actual operation of such a device. 
For instance, the design of the device ground state by 
coupling to peripheral leads, fan-out, and control of uni- 
directional logic flow are all open issues. Fabrication of 
the bowtie structure is another challenge, but recent ex- 
periments demonstrated that graphene devices only a few 
nanometers in size can be sculpted by electron beam or 
scanning probes and are quite stable 2^, 33|- Meanwhile, 
we point out that the demands on fabrication are signifi- 
cantly alleviated by the intrinsic defect-tolerability of the 
spin GNF devices: at least n — 1 out of (n -I- 2)^ — 3 car- 
bon atoms can be removed (excluding those atoms at the 
junction to which the device function is highly sensitive) 
from the dominant sublattice in an individual triangle of 
size n before the non-bonding states are eliminated and 
the local magnetic moment quenched. These consider- 
ations indicate that the devices based on the proposed 



[1 

[2; 

[4] 
[5] 

[6: 
[s: 

[9 
[10 

[11 
[12 

[13 
[14 

[15 

[16 
[17 

[18 
[19 
[20 
[21 

[22; 

[23; 

[24; 

[25; 



[26; 



p. Esquinazi et ai, Phys. Rev. Lett. 91, 227201 (2003); 
H. Ohldag et ai, Phys. Rev. Lett. 98, 187204 (2007). 
E. J. Duplock, M. Scheffler, and P. J. D. Lindan, Phys. 
Rev. Lett. 92, 225502 (2004). 

O. V. Yazyev and L. Helm, Phys. Rev. B 75, 125408 

(2007) ; O. V. Yazyev, Phys. Rev. Lett. 101, 037203 

(2008) . 

J. J. Palacios and J. Fernandez-Rossier, L. Brey, Phys. 
Rev. B 77, 195428 (2008); 

Y.-W. Son, M. L. Cohen, and S. G. Louie, Nature (Lon- 
don) 444, 347 (2006). 

O. V. Yazyev and M. I. Katsnelson, Phys. Rev. Lett. 100, 
047209 (2008). 

M. Ezawa, Phys. Rev. B 76, 245415 (2007). 
J. Fernandez-Rossier and J. J. Palacios, Phys. Rev. Lett. 
99, 177204 (2007). 

W. L. Wang, S. Meng, and E. Kaxiras, Nano Lett. 8, 241 
(2008). 

S. Bhowmick and V. B. Shenoy, J. Chem. Phys. 128, 
244717 (2008). 

A. Rajca, J. Wongsriratanakul, and S. Rajca, Science 
294, 1503 (2001). 

G. AUinson et ai, J. Chem. Soc. Perkin Trans. 1, 385 
(1995). J. Inoue et al, J. Am. Chem. Soc. 123, 12702 

(2001) . 

S. A. Wolf et al, Science 294, 1488 (2001); G. Chappert, 
A. Pert, and F. N. Van Dau, Nat. Mater. 6, 813 (2007). 
L. E. Hueso et ai, Nature (London) 445, 410 2007. N. 
Tombros et al, Nature (London) 448, 571 (2007); N. 
Tombros et al, Phys. Rev. Lett. 101, 046601 (2008). 
S. Fajtlowicz, P. E. John, and H. Sachs, Croat. Chem. 
Acta 78, 195 (2005). 

E. H. Lieb, Phys. Rev. Lett. 62, 1201 (1989). 
L. Brey, H. A. Fertig, and S. Das Sarma, Phys. Rev. Lett. 
99, 116802 (2007). 

K. Nakada et al, Phys. Rev. B 54, 17954 (1996). 
M. Wimmer et al, PRL 100, 177207 (2008) 
R. Landauer, IBM J. Res. Dev. 5, 183 (1961). 
J. M. Soler et al, J. Phys.: Condens. Matter 14, 2745 

(2002) . 

J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77, 3865 (1996). 

N. TrouUier and J. L. Martins, Phys. Rev. B 43, 1993 
(1991). 

O. V. Yazyev, W. L. Wang, S. Meng, and E. Kaxiras, 
Nano Lett. 8, 766 (2008). 

In order to obtain accurate total energies a 10~^ max- 
imum absolute change of density matrix elements was 
used as a convergence criterion for the self-consistent 
field (SCF) procedure. A Fermi-Dirac population func- 
tion (fcBT=25 meV) common for both spin-up and spin- 
down electrons was utilized. No constrains on the total 
spin of the system were imposed during the SCF runs. 
S. Pogodin and I. Agranat, J. Org. Chem. 68, 2720 

(2003) . 

[27] H. Agarwal, S. Pramanik, and S. Bandyopadhyay, New 
J. Phys. 10, 015001 (2008). 



5 



[28] C. F. Hirjibchcdin, C. P. Lutz, and A. J. Heinrich, Sci- [30] L. Tapaszto et al, Nat. Nanotechnol. 3, 397 (2008). 

ence 312, 1021 (2006). 
[29] L. A. Ponomarenko et al., Science 320, 356 (2008). 



